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Abstract 

Linear first order systems of partial differential equations of the form 
V/ = MVg, where M is a constant matrix, are studied on vector spaces 
over the fields of real and complex numbers, respectively. The Cauchy- 
Riemann equations belong to this class. We introduce a bilinear *-multi- 
plication on the solution space, which plays the role of a nonlinear super- 
position principle, that allows for algebraic construction of new solutions 
from known solutions. The gradient equations V/ = MVg constitute 
only a simple special case of a much larger class of systems of partial dif- 
ferential equations which admit a bilinear multiplication on the solution 
space, but we prove that any gradient equation has the exceptional prop- 
erty that the general analytic solution can be expressed through power 
series of certain simple solutions, with respect to the ^-multiplication. 
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1 Introduction 



We consider equations of the form 

V/ - MVg, (1) 

where /(x) and g(x) are unknown scalar-valued functions, defined on an open 
convex domain in an rt-dimensional vector space V over the field of real or 
complex numbers, and M is a constant n x n matrix. Hence, the systems under 
consideration constitute a system, that is usually overdetermined, of n first 
order linear partial differential equations (abbreviated PDEs) for two unknown 
functions. 

The Cauchy-Riemann equations 

df = dg 
dx dy 

9l = _dg 

dy dx 

are an example of an equation of the form {I}. Since there is a 1 — 1 cor- 
respondence between solutions (/, g) of the Cauchy-Riemann equations and 
holomorphic functions F = / + ig, any product 

(/ + i.9) (/ + iff) = (// - 99) + Kf9 + 9f) 

of two holomorphic functions is again holomorphic. The ordinary multiplication 
of holomorphic functions defines a bilinear ^-multiplication on the solutions 
space S of the Cauchy-Riemann equations 

*:SxS — > S 

(f,9)*(f,9) = (ff-99,f9 + 9f)- 

Moreover, since every holomorphic function is analytic, every solution of the 
Cauchy-Riemann equations can be expressed locally as a power series of a simple 
solution. For example, every solution can be formulated, in a neighborhood of 
the origin, as a power series of the simple solution (x, y) 

oo 

(f,g) = y2(a r ,b r ) * {x,y)l, where (x, y)l = (x, y) * ■ ■ ■ * (x, y), (3) 

' 7Z^— 

and real constants. 

In [B], Jodeit and Olver have given an expression for the general solution of 
the equation (TT]) . In [5] , we have have introduced a multiplication * of solutions 
for a wide class of linear first order systems of differential equations containing 
systems of the form (TTJ), where the matrix M can also have non-constant entries. 
The ^-multiplication is a nonlinear superposition principle that allows for alge- 
braic construction of new solutions from known solutions. By combining the 
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^-multiplication with the ordinary linear superposition principle, power series 
solutions can be constructed from a single simple solution. 

In this paper we show that the equations ([T]) are included in the family 
of systems equipped with a ^-multiplication on the solution space. We also 
compare power series solutions, constructed with the ^-multiplication, with the 
formulas in [5] that describe the general solution of dJ). This paper contains the 
following main results: 

1. Any equation |T]), including both the real and the complex case, admits 
a ^-multiplication which provides a nonlinear superposition formula. For 
the Cauchy-Riemann equations, this multiplication reduces to the ordi- 
nary multiplication @, obtained from the multiplication of holomorphic 
functions. 

2. Any solution of |T]) can be represented as a power series (with respect to the 
^-multiplication) of certain simple solutions. In other words, we provide 
a different and more explicit way to describe the general analytic solution 
of (TT]) . This may be compared with the Cauchy-Riemann equations where 
each solution can be obtained from the harmonically conjugated pair given 
by the real and imaginary parts of a power series in one complex variable 
©■ 

3. The equations of the form ([T|) constitute an interesting example of a large 
class of systems of linear PDEs with a simple structure, admitting *- 
multiplication of solutions. The content of this paper can be considered 
a thorough investigation of the ^-multiplication for the class of systems 
of the form H} . The results contribute to a better understanding of the 
general class of equations admitting ^-multiplication. 

Remark 1. In the special case when the matrix M consists of only one Jordan 
block for each eigenvalue, the general solution of ((7|) can be described in terms 
of component functions of functions which are differentiable over some algebra 
For this restricted class of matrices M , some of the results above can also 
be derived from the content of 

This paper is organized as follows. In section [2] we present the necessary 
results from [8] about bilinear ^-multiplication for solutions of linear systems of 
PDEs with variable and constant coefficients. Section contains a summary of 
the main results of paper [6] by Jodeit and Olver about the algebraic form of the 
general solution to the equation ([T]) . Sections Q] and [5] contain the main results 
of this paper. In section 01 we show that any equation of the form |T]) admits a 
bilinear ^-multiplication of solutions of the type described in [5]. In section [S] we 
prove (theorem [7] and theorem that every solution of |T]) can be expressed 
as a power series of simple solutions with respect to the ^multiplication. 



3 



2 Multiplication of solutions for systems of PDEs 



Certain systems of linear partial differential equations admit a bilinear operation 
on the solution space. The most general results about this operation, which we 
call ^-multiplication, are given in [5] (the results are derived for equations on a 
real differentiable manifold, but all arguments hold also for analytic equations 
over a complex vector space), and in this section we give a brief summary. 

Let Z^ = Z Q + Zip + ■■■ + Zm-in™- 1 + n m and = V + Vifi + ■ ■ ■ + 
V m -i[t m be polynomials in the variable zx, with coefficients that are smooth 
scalar- valued functions on the vector space V. Furthermore, let = Aq + 
Ax/j, + ■ ■ ■ + Aj~ii k be a polynomial where each coefficient is a n x n matrix, with 
smooth scalar-valued functions as entries. Consider then the matrix equation 

ApVVp = [mod Z M ) , (4) 

where the unknown function is a solution if the remainder of the polynomial 
A-nWp modulo Z^ is zero, i.e., if there exists a polynomial vector u M such that 

Given two solutions Vn and Wa, the *-product * Wu is defined as the 
unique remainder of the ordinary product V^W M modulo Z^, i.e., * is 
the unique polynomial of degree less than m that can be written as * W M = 
V^W/j, — QtiZ^ for some polynomial Q M . In general, neither V^W^ nor * W M 
are solutions of ©, but when A^ and Z^ are related by the equation 

A^VZ^ = (mod Z„) , (5) 

that is, when Z^ — /x m is a solution of then the ^-multiplication maps 
any two solutions into a new solution. Thus, the ^-multiplication provides a 
method for constructing, in a pure algebraic way, new solutions from already 
known solutions. Especially, when the coefficients of Z^ are not all constant, 
^-products of trivial (constant) solutions are in general non-trivial solutions. 

Example 1. If we let 

V^ = g + fit, Z^ = l + fi 2 , 

the equation Q reduces to the Cauchy-Riemann equations. Since the coeffi- 
cients in Z^ are here constant, the condition Q is trivially satisfied. Hence, 
the * -multiplication maps solutions to solutions, and it reconstructs the multi- 
plication formula (0) obtained from the multiplication of holomorphic functions: 

(/ + gn) (/ + gn) = // + if 9 + s/> + 99ft 2 

= (// - gg) + (fg + gf)ft = (/ + git) * (/ + git). 

We see that, in order to obtain the first of the main results named in the 
introduction, it is enough to find, for each matrix M, a matrix A^ and a function 
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such that the relation |[SJ is satisfied, and such that the equation ^ is 
equivalent to ((TJ) - 

The theory of ^-multiplication can also be formulated in matrix language, 
without use of congruences and the extra parameter fi. In matrix notation the 
equation ((J) can be written as 



(G) 



i=0 



where V denotes the m x n functional matrix with entries (V')ij — dVi/dxj 
and C is the companion matrix for the polynomial 



C 





1 

1 







1 



-Z 
-Zi 



o -z n . 



-Z n - 



The ^-multiplication exists for any equation ([6]) for which V = [Zq Z\ ■ ■ ■ Z m -{\ 
is a solution. The ^-product of two solutions can then be written explicitly as 
the following matrix product 



V * W = VcWcex = [zZ ViC )\zZ W t C l 



. i=0 



. i=0 



An interesting property of the ^-multiplication is that complex solutions 
can be generated, starting with simple solutions. Especially, when Z^ has non- 
constant coefficients, ^-powers of trivial (constant) solutions will in general be 
non-trivial. By using the ^-multiplication in combination with the linear super- 
position principle, one can form ^-polynomial solutions 



N 



Vn = y^a r M*; where /i* 



r=0 



fj, * fl * ■ ■ ■ * fl, 

V v ' 

r factors 



or in matrix notation 



JV 



V = ^2a r C r ei , 



(7) 



r=0 



where a, are constants. Moreover, when some conditions arc satisfied, by letting 
N — > oo in ([TJ) one obtains power series solutions. However, when the coefficients 
of Z^ are all constant (which is the case for the Cauchy-Riemann equations), 
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power series of trivial solutions will be trivial. Instead, in order to construct 
interesting solutions, one has to build power series from a non-trivial solution. 
Such power series are not considered in [5] but the results for power series of 
trivial solutions extend immediately to power series of non-trivial solutions. 
More precisely, the power series o-r(Vij) r *t which in matrix notation reads 

oo m— 1 

^a r C"ei, where C=^V t C\ 

r=0 i=0 

defines a solution of (H|) as long as the modulus of the eigenvalues of C are smaller 
than the radius of convergence of the power series ^2 a rt r an d the geometrical 
multiplicities of the eigenvalues arc constant. For power series of trivial solu- 
tions, the matrix C is identical with the companion matrix C and its eigenvalues 
coincide with the roots of the polynomial Z^. 

We can also consider *-power series with arbitrary trivial solutions as coef- 
ficients 

oo 

^2 a ^ r * O'm)*' where a^ r = a rfi + a r ,i/i H a^m-iM™ -1 - 

r=0 

Remark 2. The * -multiplication TAW!/ is a generalization of the multiplica- 
tion of cofactor pair systems introduced by Lundmark J2J/ ; which in turn is a 
generalization of a recursive formula for obtaining new cofactor pair systems 
I VA \10f . A cofactor pair system is a Newton equation q h + T^cfqi = F , 
h = 1, 2, . . . ,n on a Riemannian manifold, where the vector field F can be writ- 
ten as 

F = -(det J)~ 1 JW = -(det J)" 1 JVV 

for some functions V, V and two special conformal Killing tensors J , J. The 
family of cofactor pair systems contain all separable conservative Lagrangian 
systems, that in general can be integrated through separation of variables in the 
Hamilton- J 'acobi sense ]1'A Ub\ IT71 [751 [3J [JJ/. The multiplication of cofactor 
pair systems is a mapping that, for fixed special conformal Killing tensors J 
and J, for two solutions (V, V) and (W, W) of the equation (det J) _1 JVF = 
(det J)~ 1 JW prescribes a new solution (V, V) * (W, W) in a bilinear way. 

3 The general solution of V/ = MVg 

In [B], Jodeit and Olver describe a general analytic solution of the equation (I} 
for any constant, real or complex, matrix M . In this section, we introduce some 
notation, and briefly describe the main results from [5J. 

A linear change of variables x — * Ax transforms the equation ([T|) into V/ = 
A~ T MA T Vg. Thus, by performing this kind of transformation, we can assume 
that the matrix M is in some suitable canonical form with respect to similarity. 
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We will assume in the following that M is in Jordan canonical form when V is 
a vector space over the complex numbers, and in real Jordan canonical form [5] 
for real V. 

3.1 The complex case 

If we let Ai, . . . , X p denote the distinct eigenvalues of the matrix M, there is a 
primary decomposition of the vector space 

V = V 1 © • • • © V p (8) 

consisting of invariant subspaces V k = kcr (M — Xkl) ek , where efc is the geo- 
metric multiplicity of Afc . As mentioned above, there is no loss of generality to 
assume that M has a diagonal block structure M = diag(A/ 1 , . . . , M p ), where 
M = diag(J^, . . . , J k k ) and each J? is a Jordan block corresponding to the 
eigenvalue Afc, i.e., 

Afc 1 

Afc 

'•■ 1 

Afc 

There is a simple decomposition of the general solution of ([1]) that is associated 
with the primary decomposition ([5]) of the vector space V. Every solution (/, g) 
of the equation JT]) can be written as a sum 

f = f+f' 2 + --- + F, 9 //' • <f 9 P , (9) 

where (f k , g k ) is a solution of the corresponding equation V/ = M k Vg. Hence, 
the problem of describing the general solution of ([I]) is reduced to the case when 
the matrix M consists of Jordan blocks corresponding to a single eigenvalue 
A. Therefore, let M = diag(Ji, . . . , J m ), where Jfc is a Jordan block of size 
(nfc + 1) x (rifc + 1) corresponding to A. We assume also that coordinates are 
chosen in such a way that the size of the Jordan blocks is decreasing, i.e., n\ > 
ri2 > • • ■ > n m > 0. Any vector x G V is decomposed as x = [x 1 x 2 ■ ■ • x m ] T 
so that x k = [xq x\ ■ ■ ■ x k lk ] T contains the variables that correspond to the 
block Jfc. The variables x\, Xq, . . . , X™, that appear at the first position in the 
vectors x 1 , x 2 , . . . , x m , are called major variables and all other variables are 
referred to as minor variables. For each block Jfc, we introduce a function 

x k =x k ) +x k lf i+--- + x k k fi nk , 

which depends on the variables contained in x fc and a real parameter fj,. Fur- 
thermore, let Vk denote the number of Jordan blocks of M which are of size at 
least k + 1, i.e., — max{i S N : rij > fc}, and define vectors x|j , . . . , x^" 1 ^ in 
the following way: 

x (fe) _ r r l r 2 ... T " k ] T 
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According to [B], two analytic functions /(x) and g(x) form a solution of the 
equation |T]) (where we have assumed that the matrix M has only one eigenvalue 
A and has a Jordan canonical form) if and only if 

/ = fx + /2 H h + c, 5 = .91 + g 2 H h ff ni , 

where c is an arbitrary constant, and there exist scalar- valued analytic functions 

4>k(s 1 ,...,s Vk ), k = 0, 1, . . . ,m, 

such that /fe, ^fe are given by 



5H X 



<9 fe 

Qk 

dfi k 



AW = A^fc (x« 
& (*«) 
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k-l 



k-l 



M=0 



(10) 



There are two interesting observations regarding the general solution: 

1. By a change of dependent variable h = f — Xg the equation |T]) transforms 
into an equation for h and g which is independent of the eigenvalue A. 



The functions (f>o,(j)i,... ,4> ni in pop that are arbitrary, depend only on 
the major variables Xq, Xq, . . . , Xq 1 , whereas the dependence on the minor 
variables of the functions and g^ is restricted to certain fixed polynomi- 
als. For example, in the case when m = 1, i.e., when the matrix consists 
of only one Jordan block, the functions fk and gu can be written as 



k-l 

fk{x) = Agfc(x) + k} j Pk-i,j(xi,x 2 , 

.,x k - j+ i)<f) k 3 \xo), 



x k-j)(t>k\x ) 



9k 



x 2 , 



(i) 

where (jr k denote the j-th derivative of the function and Pk.j are fixed 
polynomials, related to the partial Bell polynomials [2]. 

We provide a concrete example in order to illustrate how to read this result 
from [Hj. 

Example 2. Consider the equation ([7]) where M is a constant 5x5 matrix. 
As mentioned above, by changing independent variables in a linear way, we can 
assume that M is in canonical Jordan form. We consider the equation ([7]) in 
the particular case when M is given by 



M 



Ai 1 

Ai 











A 2 








1 

A 2 









A 2 
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where Ai and A 2 are distinct eigenvalues. First of all, the general solution can 
be decomposed as f = f l + / 2 , g = g 1 + g 2 , where (/ , g 1 ) and (/ 2 , g 2 ) are 
solutions of 



Ai 1 
Ai 



V5 1 , and Vf 2 = 



A 2 1 
A 2 



A, 



and 



respectively. From the formula we then obtain 

\ 9 1 = 9o +9i= Xi4>'i(xo) + 4>o(x ) ' 

\ g 2 = 9o + g\ = to(yl, 2/0) + yWiivl) ' 

where xo, Hq, Uq are the major variables, x\, y\ the minor variables, and 
4>q, 4>i, ipo, ipi are arbitrary scalar-valued analytic functions. 



3.2 The real case 

We consider now the equation |T]) in a convex domain of a vector space over 
the real numbers. Since the field of real numbers is not algebraically closed, the 
matrix M is in general not similar to a matrix of Jordan canonical form. Instead 
we assume that M — diag(M 1 , . . . , M p ) has the real Jordan canonical form, so 
that the blocks corresponding to the real eigenvalues have the same structure as 
in the complex case, whereas a block M k — diag(Af^, . . . , M* ) corresponding 
to a complex conjugated pair of eigenvalues, A, A = a ± consists of blocks 
Mj^ of the form 

" A I 

A '■■ 

'■■ I 

A 

The general solution of |1| can again be decomposed as in ^ , where each pair 
{f k >g k ) is a solution of the subsystem V/ = M k Vg. Thus, since the equations 
corresponding to real eigenvalues can be treated in the same way as in the 
complex case, it is enough to consider the case when M = diag(Mi, . . . ,M m ) 
has only one pair of complex conjugate eigenvalues A, A. We assume that each 
Mi is a matrix of the form (fTTj) of size 2(m + 1) where n\ > ■ ■ ■ > n m > 0. 
Moreover, we denote the variables corresponding to the Jordan decomposition 
by £0 1 Vo > ■ ■ ■ > x ln >Vni ' ) vl > ■ ■ ■ > x n m , Vn m > and denne complex variables = 



, where A = 



a (3 
-f3 a 



I = 



1 
1 



(11) 
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x^+iy!?, z'f 




ii/j. In a similar way as in the complex case, we introduce 



functions 



z 



,fc 

7< 



A; = 



1,2, ... ,771, 



depending on a real parameter /x, and vectors 




2^ fc ] , k — 0, 1, ... , 



where 7^ is defined in the same way as above. If we let F = f — Xg, the general 
analytic solution of V/ = MVg can be expressed as F = Fq + F\ + ■ ■ ■ + F ni 
where 



and <f>k(si, ■ ■ ■ , s Vk ) are again arbitrary complex- valued analytic functions de- 
pending on Vk complex variables. 

3.3 A characterization of the general solution in terms of 
different iable functions on algebras 

In [14] , the general solution of ([!]) is characterized through components of func- 
tions which are differentiable over algebras. From this characterization, some 
results regarding multiplication of solutions can be obtained. However, the re- 
sults in P3] , regarding the equation (fTJ) is given only for a quite restricted class 
of matrices M . We will give a brief summary of the concepts and results in [14] 
about the equation ([1]). 

Let A be a commutative algebra of finite dimension over the real or complex 
numbers, and consider a function / : A — * A. We say that / is A- differentiable 
in the point a G A if the limit 



exists. When A is a C-algebra, any A-differentiable function is analytic, i.e., 
it can be expressed locally as a power series in the A-variable x. A function 
V : A — * R (or C) is called a component function of an A-differcntiable function 
/ if it can be written as V = <fi ° f for some linear function <fr ■ A — > R (or C). 

When a 77, x 77, matrix M consists of one single Jordan block, there exists 
an algebra A which is generated by a single element s such that in some basis 
&i, . . . , a„, M is the matrix of multiplication by s, saj = ^ M^aj. If V = <fi(f) 
is a component function of a C 2 A-diffcrcntiable function / and W = <f>(sf) is 
the corresponding component function of sf, then \7W = MW . Conversely, 
if V and W are C 2 functions such that VW = MW, then V = <j>{f) and 
W — <fi( s f) f° r some A-diffcrcntiable function / and some linear function <f>. 

This characterization in terms of component functions of A-differentiable 
functions implies the existence of a multiplication of solutions for the equation 




(12) 



hm x- 1 (/(a + x)-/(a)) 
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|T]) since products of j4-differentiable functions are again A-differentiable. More- 
over, in the complex case, since any A-differentiable function is analytic, one 
can prove that any analytic solution of (TTJ) can be expressed as a power series 
of a simple solution. Thus, since some of the problems we are considering in 
this paper are already treated in [14] , it is worth to emphasize in which way we 
extend these results. 

1. In [14j . the only case considered is when the matrix M has exactly one 
Jordan block corresponding to each eigenvalue. In this paper, the equation 
(H|) is studied for all matrices M. 

2. In [T3], only the existence of the algebra A, from the characterization of 
the general solution of (JTJ) , is given. An explicit construction is natural 
and simple in the complex case, but not trivial in the real case. In this 
paper we provide an explicit multiplication formula for solutions of (TIJ for 
every possible M , both in the real and in the complex case. 

3. The characterization of ^4-differentiable functions in terms of analytic func- 
tions is only valid for the complex case. We show that any solution of ((T|), 
in both the complex and real case, can be expressed as a *-power series. 

4. The equations {1} constitute only a simple special case of a quite large 
family of systems of PDEs with ^-multiplication. Therefore, the results 
obtained in this paper, for the restricted class of systems |TJ, indicate 
that similar results may be established for the entire class of systems with 
^-multiplication. 

In this paper, we have used the notation and conventions used in [6] and 
[8], rather than in [14]. However, for the special cases when a multiplication of 
solutions can be obtained from |14j . this multiplication coincides, even though 
it may not be obvious immediately, with the ^-multiplication studied in this 
paper. 

4 Multiplication for systems V/ = MVg 

In this section we show that every equation |T]) admits a ^-multiplication. The 
approach we use for proving this, is to extend the equation ([1]) by introducing 
certain auxiliary dependent variables and adding some equations that are conse- 
quences of the original equation. One can then show that the extended system 
has the form ([4| and that it satisfies {5]), so that it admits a ^-multiplication. 
Since there is a simple correspondence between solutions of the original and 
the extended system, one can interpret the ^-multiplication as an algebraically 
defined operation on the solution space of the equation |T]). We treat the real 
and complex case separately, starting with the latter. 

Remark 3. We will sometimes use the exterior differential operator d instead 
of the gradient operator V. Since we are only considering complex analytic 
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functions, the exterior differential reduces to the d-operator J^. Thus, we can 
write the equation ([7]) as df = M T dg which should be understood as the following 
equation for row matrices: 

[dif • ■ ■ d„f] = [d l9 ■ ■ ■ d n g]M T , 

where di = d/dxi denotes the partial derivative with respect to Xi, and M T 
denotes the transpose of the matrix M. 



4.1 The complex case 

When the general solution of the equation (fT]) is obtained in [5] , Jodeit and Olvcr 
split their investigation into three cases, and we will use the same strategy in 
this paper. The cases are: 

1. The matrix M consists of a single Jordan block, i.e., 

A 1 



M = 



A 



1 
A 



(13) 



2. M = diag(Ji, . . . , J m ) has only one eigenvalue but consists of several 
Jordan blocks corresponding to that eigenvalue. 

3. M is a general complex matrix. 

Obviously, the first case is a special case of the second case, which in turn is a 
special case of the third case. However, treating the first two cases separately, 
has a didactic advantage. 



4.1.1 One Jordan block 



We assume now that the matrix M consists of one Jordan block of size n + 1. 
By a change of dependent variable h = f — Xg, the equation (JTJ) reduces to the 
equation 



V/i = U n Vg, 



(14) 



where U„ is defined as the the Jordan block of size (n+1) X (n+1) corresponding 
to the eigenvalue A = 0, i.e., 



U n := 



1 



1 




(15) 



Equation (|14p admits a ^-multiplication on the solution space: 
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Proposition 1. Let B be any nilpotent constant matrix. Then the system 
V/ = BSJg can be extended to the system 

A M VV^ = {mod fi r ) , (16) 

where 

A^ = -B + nl, Vp = V Q +nV l +>-- + //" 3 v;-3 + S'^f + ^9, 

and r is a natural number such that B r = 0. The functions Vo, ■ ■ ■ , V r s are 
auxiliary functions that are uniquely determined (up to additive constants) by 
any given solution (/, g) of \H\ . 

The proof relies on the following lemma. 

Lemma 2. Let B be a constant n x n matrix, and f any smooth function such 
that d(Bdf) = 0. Then d(B k df) = for any keN. 

Proof. The proof is by induction over k. For k — 2, by the assumption d(Bdf) — 
0, we have 

n 

(d(B 2 d/)) y . - {BabBbidadjf -B ab B bo d a d % f) 

a, 6=1 

71 

= (B aj B bi d a d b f - B ai B bj d a d b f) = 0, 

a, 6=1 

where Bij denotes the entry in row i and column j of the matrix B. Assume now 
that the statement is true for k = p—1, where p > 3. If d(Bdf) = 0, there exists, 
according to the Poincare lemma, a function g such that dg = Bdf. Thus, we 
have d(Bdg) = d(B 2 df) = 0, and therefore by the inductive assumption 

= dtB^ds) = d(B p df). 

□ 

Proof (of proposition^. Let — —B+fil, and — Vo+/iVi + - • •+fi r ~ 1 V r —i. 
Then, the components of equation (|16|) read 

= B T dV 
dV a = B T dV 1 

< 

dV r ^ 2 = B T dV r - 1) 
or equivalently (using the assumption that B is nilpotent), 

dVi = (fl r ) r-1 ~ < dV r _i, i = 0,l,...,r-2. (17) 
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It is then obvious that if solves (|T71) . then the functions / = V r -2, g — V r -\ 
solve the equation V/ = BVg. 

On the other hand, assume now that (/, g) solve the equation V/ = BVg, 
and let V r - 2 = /, V r -i = g. Then, since d(B T dK_i) = d 2 V r - 2 = 0, the 
previous lemma implies that d((B T ) k dV r -i) = 0, for any k > 1. Therefore, 
according to the Poincare lemma there exist functions Vq, Vx, . . . , V r -3 such 
that (fTTI) is satisfied. □ 



Since A^Vfi 7 ' = is trivially satisfied, the equation (flT)]) admits a ♦-multi- 
plication. Given any two solutions and V^, of (|16[) . their ^-product is given 
by the formula 



fc=0 \j=0 



or using the matrix notation 



^0 












Vi 












. V r -i _ 




. V r -l _ 





v v 



VoV r -l + ViVr-2 



+ V r -lVo _ 



(18) 



Especially, since the matrix U n is nilpotent with £/™ +1 = 0, proposition [T] allows 
us to extend the system (fT4|) to the ^-dependent system 

(19) 



{-U n + /i/)VV r M = (mod n n+1 ) 



Thus, given two solutions (h,g) and (h, g) of (fT4|) . there exist unique (up to 
additive constants) functions Vq, . . .V n and Vq, . . . V n such that 



W i = UZ- i Vg, and VV i = U^- i Vg, 



0,1, 



A new solution of (fT4")) is then algebraically constructed through the ♦-multi- 
plication as 



(n-\ 



(Kg) = $>*U-<, Z!^"- 4 



\i=0 i=0 / 

4.1.2 Several Jordan blocks corresponding to one eigenvalue 

Assume now that M = diag(Ji, . . . , J m ) has only one eigenvalue A, where each 
Ji is a Jordan block of size rij + 1, and ni > ri2 > • ■ ■ > Just as in the case 
of one single Jordan block, by changing one dependent variable h = f — Xg, we 
obtain the equivalent equation 



V/i = UVg, where U = diag(L/" ni , . . . , U nm ), 



(20) 
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and Ui is given by (fT5|) . Since U is nilpotent, we can apply proposition [1] with 
B = U and r = n\ + 1, in order to extend ([20]) to the system 



(-U + /LtJ)V^ ee (mod , (21) 

Thus, there is a ^-multiplication for solutions of the equation ([T]) (having the 
same multiplication formula (fl8|) ^ also for the case when M consists of several 
Jordan blocks corresponding to the same eigenvalue. 

Besides generating solutions of the system (|20p with ^-multiplication, it fol- 
lows from the next proposition that one can also embed solutions of the subsys- 
tems V/i = diag(J7„ ; , . . . , Un m )Vg, for any 1 < i < m, into the solution space 
of ([20]). 

Proposition 3. Suppose that is a solution of the equation 

r-i 

A m W m eeO (mod//), V li =^V i fj, i ( 22 ) 

defined in a subset Cli of a vector space V\, where = —A + fil, and r is a 
natural number. Then — fJ. s V^ is a solution of the extended system 

r+s-l 

B m VW m eeO (mod //•+*), W ^ 

defined in Qi x V2, where V2 is another vector space, = B + fil , and B is 
the square matrix with the block structure 



B 



C 
A 



for some constant matrix C and some s £ N. 

Proof. Let be a solution of the equation (|2"2"|) , i.e., there exists a 1-form a, 
not depending on /i, such that A^dV^ = fi r a. Then 



Bjd( M %) = //SjdV; = M s (CjdiV; + Ald 2 Vj 



AtM^d 2 V M = M s+r a = (mod y r+s ^ 



where di and d 2 denote the exterior differential operators on V\ and V2, respec- 
tively. □ 

Let (h l ,g l ) be a solution of V/i l = diag(C/ ni , . . . , U nm )Vg l . Then, according 
to proposition [U there exists a unique (up to unessential constants) solution 
V* = Vtf + h V* i _2lt ni ~ 2 + h 1 ^" 1 + g l n n > of the extended system 

(-diag(C/ i; ...,U m ) + fxl) dV; ee (mod ^ +1 ) . 

According to proposition^ it then follows that V* can be embedded as a solution 
H V' of the equation ((20]). 
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4.1.3 The general complex case 

Since the general solution of equation (Q]) can be written as a sum © where 
{f l ,g l ) is obtained from the general solution of the system ([20]) , it is clear 
that there is no non-trivial ^-multiplication for the system |T]) other than the 
multiplications induced from the subsystems corresponding to single eigenvalues, 
that was introduced above. 

In order to illustrate the ^-multiplication for equations of the form ([T]), we 
return to example [5] 

Example 3. Consider again equation ([7]) in the case when M is the 5x5 matrix 
defined in example^ Take for instance the particular solution 

f = \x2xqXx + X^xof + (x ) 2 + Xiylivl) 2 + A 2 2/i + Va 
g = 2xa Xl + {xq) 3 + yKvlf + y\ 

This solution can be written as a sum of * -products of simple (first order poly- 
nomials) solutions of subsystems. In detail, 

f = Xigi + hi + A 2 (?2 + h 2 
g = 9i + 52 

where 

hi + mn = l l * (xa + l iX i)l + ( X Q + V-Xi)l, 

and 

h 2 + m = (KvDl) * ivl + wi) + (vo + wl)- 

We note that the factor ^(j/g)^ is the embedding, according to proposition^ 
of the solution (jJq) 2 , — (Ho) 2 of the subsystem Vh — U$Vg (or equivalently 
h'iVo) — Q)> into the system \/h 2 = diag(J7i, Uo)Vg2- 

The possibility to express solutions of an equation of the kind as a sum 
of simple *-products, that was illustrated in the previous example, is a gen- 
eral property. In the next section we show that any analytic solution of |T]) 
can be expressed as a power series of simple solutions with respect to the *- 
multiplication. 

4.2 The real case 

We consider now the equation |T]) over a real vector space when the matrix 
M has real entries. From the discussion in section 13.21 it is clear that it is no 
restriction to assume that M has only one pair of complex conjugate eigenvalues. 
Thus, we let M — diag(Mi, . . . , M m ) be an arbitrary square matrix having real 
Jordan form with eigenvalues A, A = a ± i/3, and we use the same notation as 
in section 13.21 Furthermore, we can, without loss of generality, assume that 
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a = and /3 = 1 . This is realized by changing both dependent and independent 
variables according to 

/ = / - otg, g = (3g, 
~x) = ^x), yk =/3 3-n ky k^ i = ,l,...,n fe fc = l,2,...,m. 

The matrix M is not nilpotent, which was the case for the complex case, but 
the system (p} is a so called quasi-Cauchy-Riemann equation which are known 
to admit ^-multiplication. 

Theorem 4. Let M be a square matrix of size 2(n + l) having real Jordan form, 
and with no other eigenvalues than ±i. Then equation (fl|) can be extended to 
the system 

A M W M = (mod (fx 2 + 1)™ +1 ) , (23) 

where 

= M- 1 +fil, Vfj, = Vq + Vifi + • • • + V 2n+1 ^ n+ \ 

and Vq = f and V^n+i — 5- The functions V\, . . . , Vi n are uniquely determined 
(up to additive constants) by the solution (f,g)- 

Proof. Since detM = (det A)™ +1 = 1, equation (fTJ) can be written as 

M- 1 V/ = detM- 1 V 5 . (24) 

Equation (|24l) is a quasi-Cauchy-Riemann equation and is therefore equivalent 
to the parameter dependent equation (see [7] J]]: 

(AT 1 + fiI)VVf, = (mod det (M _1 + fj,I)) , (25) 

where Vq — f and Pssn+i = 9- Thus, since 

det (M- 1 + fil) = (dot (A + nl)) n+1 = (/i 2 + 1)" +1 , 

the proof is complete. □ 

Since (M _1 + /i/)V(^ 2 + l) ,l+1 = 0, the solution space of the system (g3J 
admits ^-multiplication. In the special case when (fTJ) is the Cauchy-Riemann 
equations, i.e., M = A with a — and [3=1, the parameter dependent equation 
(|25"|) coincides with the original equation {1} and the ^-multiplication reduces 
to the ordinary multiplication, obtained from the multiplication of holomorphic 
functions. 

In a similar way as in the complex case, solutions of (fl]) can be constructed 
by embedding of solutions of subsystems. Let V^ be a solution of 

i^VV^ = (mod ( M 2 + 1)" +1 ) , 

where A^ = diag(M i ~ , . . . , M^ 1 ) + fil and h = m + ■ ■ ■ + n m +m — i. Then 
(/i 2 + l) n ~ n V^ is a solution of ([23)1 . A proof of a more general result is given in 

i- 

1 In [7], the definition of a quasi-Cauchy-Riemann equation require that the matrix M is 
symmetric, but as long as the entries are constant this assumption is not necessary. 
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5 Analytic solutions are *-analytic 



We have seen in the previous section that any equation of the form (QJ can be 
reduced to a number of systems of the form (f2"Tj) (or (|23[) in the real case) , each 
one equipped with a ^-multiplication on the solution space. In this section we 
prove that every analytic solution can be expressed as a power series of simple 
solutions, similarly as for the Cauchy-Riemann equations. Being precise, by 
simple solution we mean a solution which is linear in the independent variables. 
Furthermore, we say that a solution (/, g) of |T]) is ^-analytic if it can be ex- 
pressed locally as a finite sum of power series of simple solutions with respect to 
the different ^-multiplications of the corresponding subsystems of the form (|2ip 
and (|23p . In other words, in this section we aim to prove that every analytic 
solution of ([T]) is also ^-analytic. Since it is no restriction, we will only consider 
power series expansions in a neighborhood of the origin. 

In both the real and the complex case, equation |T]) is extended to a finite 
number of /x-dependent equations of the form 

(X + M/) W M = (mod Z„) , (26) 

where X is a constant matrix, = Zq + • • • + Z n fi n is a polynomial with 

constant coefficients, = Vq H 1- V n fj, n , and V n = g. Equation (|26]) can be 

written without the parameter /x as a system 

XVVj + Wj_i = ZjVg j = 0,l,...,n, (27) 

where V-i = 0. Let V/[ N ^ = Vq + • • • + Vn N ^ fj, n be a solution with coefficients 
being polynomials of degree N in the coordinate functions, 

JV 

fc=o j=(to,...,t n )eN n 

ioH \-i n = k 

Because of the relation (|2"T|) . the power series 

lim V 3 (N \ j = 0,l,...,n 

will all have the same domain of convergence. This means that when construct- 
ing *-power series of a simple solution, it is enough to consider the domain of 
convergence for the power series defined by the coefficient at the highest power 



5.1 The complex case 

We split the investigation of the complex case into three different cases in the 
same way as in the previous section. 
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5.1.1 One Jordan block 

We consider first the equation ([l} when the Jordan decomposition of the con- 
stant matrix M consist of one single Jordan block. We have seen that, then 
there is no loss of generality to consider the equation (fH|) instead. 



Theorem 5. Let (h,g) be an analytic solution of the equation Then there 

exist constant solutions (bi)^ — & 4 o + bn/i + • ■ • + bi n /i n (with bij £ C) of the 
extended system (j 19h for i = 0, 1, . . ., such that 

oo 

* Ml = V + ■ ■ ■ + K-2// 1 - 2 + h^- 1 + 9f i n , (28) 

2 = 

where = xq + Xifi + • • • + x n \x n . 

We formulate a result about the structure of the *-powers {x^)\ as a lemma: 
Lemma 6. For each i G N we have 

n 

M=5>. M ( X )^, (29) 
j=o 

where are the polynomials given by 



*M*) := E 



ao+oiH ha„=i 

ai+2a2 4 \-na n —j 



a®, a\, . . . , a„ 



ao„ai . . . „a„ 
x •''l x n 



= E ( ) X "-SjsC^l) 2a; 2, ■ ■ ■ , {j-S + iy.Xj-^i), 

where Bj s = Bj S {x\, Xi,..., Xj-^-i) are the partial Bell polynomials Jl|/ given by 

,. ,3 . -^A^ w •A(i- s+ i) ! ; 

an+aiH haj-iH-l=s 

ai+2a 2 H h(j-»f l)a3-»|-l=j 

Proof. From the multinomial theorem, we have 

4= E 



ao,ai,...,aTi€N 
ao+axH ha„ — i 



i 

ao, flti, . . . , a, 



„, a ™ a i ™Ora ,,ai+2a 2 H kia„ 

•''O x l ,J, n r 



Thus, formula (|29|) follows now from the definition of the ^-multiplication. More- 
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over, 



•«w - E |r E 



a =0 



E 



ai ,02 , . ..,a n (EN 

diH hcin.— 2— do 

ai+2a2H h^a T i — J 



ai!a 2 !---a n ! 1 



ti 



E 



ai,a2,...,aj_ £ 4.iGN 

aiH haj-s+l = s 

a 1 +2 a 2 H h ( J-s-f-1 ) a j-s+i = J 



«1, 



E ( ) x o "-Bia^i^ 2 ^, • • • , (j-s + iy.xj-^). 

s=0 W 



□ 



Proof (of theorem^). For any solution (h,g) of (|14jl . /i is uniquely determined 
from g (up to an additive constant). Thus, it is enough to prove that for 
any analytic solution {h,g), we can choose {bij^ such that the highest order 
coefficient of X^o^iV * Ov)* coincides with g. From the results obtained in 
[6], presented in section [3J the general solution g is given by 



9( x ) = E jrj fa (*/*) 

3=0 °^ 



(30) 



M=0 



where 0i, <^>2> • • • > 4>n are arbitrary analytic functions of one complex variable, 
given by 



i(s) =^2,c%js\ j = 0,1,..., 5 



(31) 



Hence, if we write the function g in a more explicit way, we obtain 



j=0 i=0 V ' 



EE' 

j=0 i=0 



duJ 

^ a ,...,a„SN 
\ a + ---+a 7l = 



E 



M=0 



ao, . . . , a, 



l x 1 



EE 

j=0 i=0 

We define now the constant solutions (60 m as 



T «n ,,oiH kia„ 



= E( n - fc ) 
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From lemma §§§ , it then follows that the highest order coefficient of the *-power 
series — Y^(bi)n * ( x n)* coincides with g(x). 

□ 

We illustrate this result with a simple example. 

Example 4. Consider the system Vft. = C/2V 'g, where U2 is the 3x3 matrix 
described above, i.e., 



V/i = 



1 
1 




v 9 . 



(32) 



The general analytic solution of \32\ is given by 



h = <Pi(xq) + 2xi(f>' 2 (x ) + c 

g = c/)o(xq) + xi4{(xo) + 22:2^2(^0) + xf<f>2(xo) 



(33) 



where 4>o, 4>i, 4>2 are general analytic functions of one complex variable. If the 
corresponding power series representations are given by (| <?-?[) for j = 0, 1, 2, then 
the general solution (jffgp is ^-analytic with the following power series represen- 
tation: 



' V ' 


OO 


Ci2 




x 


h 




Cil 




X\ 


9 


i=0 






. X2 



5.1.2 Several Jordan blocks corresponding to one eigenvalue 

Assume now that the matrix M in equation {T]) consists of only one eigenvalue 
A, but that its Jordan canonical form consists of several Jordan blocks. As we 
have seen, it is then no restriction to assume that the eigenvalue is zero, so 
that (fT]) is reduced to the equation (|2"0"|) . Then a given analytic solution can in 
general not be written directly as a power series of simple solutions with respect 
to the ^-multiplication (TTg|) for that system. This is easily realized from the 
following example. 



Example 5. Consider the equation Vh = diag(£7i, J7o)V<7, i.e., 

"010" 

Vh= Vg. 


The general analytic solution is given by 



(34) 



h = (j>(x) + c 

9 = y<f>'(x) +i>(x,z) 
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where <f> andip are arbitrary analytic functions, and the * -multiplication formula 
is given by 



(35) 



" h ' 




' h ' 




hh 


9 


* 


9 




Kg + gh 



Thus, we see that the ^-product of polynomial solutions will never have higher 
degree in the z-variable than the highest degree of g and g. Since a polynomial 
solution can have any given polynomial degree in z, it becomes obvious that 
not every analytic solution of \34\ can be written as a power series of simple 
solutions with respect to the * -multiplication <\35\ . 

The remedy is to embed power series solutions of the trivial subsystem V/i = 
Uq Vg in the sense of proposition [3J The general solution is given by g = t(z) 
where r is an arbitrary function, and the ^-multiplication coincides with the 
ordinary multiplication. Thus, according to proposition [3J we can embed it- 
powers zi — z° of the trivial solution g — z into solutions of the system \34\ - 
Any analytic solution of <\34\ can then be written as an infinite sum of products 
of the solutions (x, y)\ and the embedded solutions (0, z J ) of the ^-powers zi: 



i=0 

oo oc 

^iciyx 1 + 2J (Hf 

i,3=0 



.* ?3 



i=0 



i=0 



i=0 



In the general case with one eigenvalue, we proceed in the same way as 
in the example, and embed power series solutions for the subsystems V/i* = 
diag(J7 ni , . . . , Un m )X7g l , in the sense of proposition |3J 

Theorem 7. Every analytic solution of the matrix equation <\2U§ is ^-analytic. 
In detail, let (h, g) be an analytic solution defined by 

m 



E 



3=0 



Li) 



fi=0 



where (f>o, 4>\, ■ ■ ■ , 4> ni are analytic functions with power series representations 

<f>j(ai,32,...,8 Vj ) = 



E 



Cj f <j 5-^52 * * * 



Then (h,g) is ^-analytic with power series representation 

m 

Vb + --- + K 1 -2M" 1 - 2 + ?iM ni - 1 +5M" 1 =E E 

r=l Z=(i I ,i2,...,i r )6N r 

ir), is the trivial solution of 



where, for each multi-index I = (ii,«2>- 
the extended system <\21l defined as 

j=n T+l + l 



(n 



m+1 
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(*„)£ == * (M ni " ri2 (4)? * • • • * (^-^W)) . (36) 

Remark 4. TTie definition <\S6} of the powers (x^)^ needs some further explana- 
tion since the ^-operators on the right hand side of \3€ty denote in fact multiplica- 
tions for different systems. In order to calculate (x^)£, one has to start with the 
power (xV) % J , which is a solution of the subsystem V/i = diag([/ r , . . . , t/ m )Vg. 
This solution is then embedded, in the sense of proposition^ into the solution 
pTh—i— nr {x r y^ of the subsystem V/i = diag([/ r _i, U r , . . . , U m )Vg, and is there- 
after ^-multiplied with (a;^ -1 )/ -1 . One continues this successive embedding of 
solutions into larger subsystems until the solution (a; 2 )* 2 * • • • * (/i" r_1 ~™ r (a^)i r ) 
of the system V/i = diag(t/2, ■ ■ ■ , U rn )Vg is embedded (through multiplication 
with fi 7ll ~ n2 ) into a solution of the system and then finally multiplied with 

the solution (x*)* 1 . 

As in the case with only one Jordan block, we place the calculation of the 
powers (x M )^ in a separate lemma. 

Lemma 8. For each I = (ii , 22, ■ ■ ■ , v) S N r we have 

n r 
3=0 

where are polynomials defined through the polynomials ^i.j in lemma\^ as 

* 7j (x 1 ,x 2 ,...,x'-)= ^.i! (x 1 )*i 2)ia (x a )-..* Wr (x'-) 

0<ji,...,ir<n r 
ji+iaH hir=j 

Proof. For r = 1, the statement reduces to lemma |6l For arbitrary r, since 
. . , x r ] T , we may assume by induction that 



(x^ = £* il;S (xV * U ni -" 2 ^% 2 ,... ir ), t (x 2 ,...,x'-K 



,n„*>2-T«r+t 

Vs=o / V t=o / 

Thus, from the ^multiplication formula (|18[) . we obtain 

( X ^E (E^ 1 ^( xl )*fe,...^) J -,(x 2 ,-..,x'')) 

n r 

= £* JJ (x 1 ,x 2 ,...,x>^-«^, 

where the last equality follows immediately from the definition of the polyno- 
mials i&ij. □ 
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Proof (of theorem^. Since v nr+1 +i = ^n r+ i+2 = ••• = v Ur = r for r — 
1,2,..., to, we have 



= £ £ E^^j^^-Kr) 



Thus, it is enough to show that for each J = (ij., 12, . . . , v), the expression 



E ^a7j(W 1 (^) i2 ---K)") 

j=n r+1 +l 



coincides with the highest order coefficient of * that, according to 

the previous lemma and the ^-multiplication formula (fl8|) . equals 

E j' !c ^*^- 

j'=n r+ i + l 

The equality is therefore proven by the following calculation: 



^((»ir(«Sr-(*y) 



/x=0 



s=l 

imi-{ Vi r n T 



M=0 



= h E W 



fj,=0 



□ 



5.1.3 The general complex case 

Since the general solution of |T]) can be decomposed as a sum © of solutions to 
subsystems, one of the main results of this paper is an immediate consequence 
of theorem 

Theorem 9. Over the complex numbers, every analytic solution of the matrix 
equation (fl|) is *- analytic. 



5.2 The real case 

The sub-equations of ([1]) corresponding to real eigenvalues of M are treated in 
the exact same way as in the complex case. Therefore, it is enough to consider 
equations corresponding to complex conjugate pairs of eigenvalues. 
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We consider only the case when M consists of a single real Jordan block, 
i.e., M has the form (TTTj) . As we have seen in section[4j there is no restriction to 
assume that a = and (3 = 1. We have already proved that the corresponding 
system ([T]) can be extended to the parameter dependent equation (f2"3")l . 

In the complex case with one single Jordan block, there was a very natural 
choice of "simple" solution to construct *-power series form. Namely, x M = 
xa+xi[i+- ■ -+x n \i n where the coefficients are given by the coordinate functions. 
In the current situation on the other hand, the corresponding function, with 
coefficients given by the coordinate functions, is not even a solution. Instead 
of searching for an appropriate simple solution, suitable for constructing power 
series, we change variables x = Bs, where B is a constant matrix, in such a way 
that the simple function 



so 



2n+l 



sifi-\ h s 2 «+i/i 

becomes a solution. 

Lemma 10. s M is a solution if B is given by the block structure 



(37) 



£ = 1,2, 
3 = 1,2, 



.,n + l 
.,2(n + l) 



1 " 




' 1 " 


-1 


, ei = 






where Bij denotes the 2x1 block in the block-position of B, and 
A = 

Proof. By changing variables x = Bs, equation (|23[) transforms into 
(B T Ar 1 B- T + fiZ)W M = (mod (/i 2 + l)' i+1 ) . 

Z 



If we let X = B T M- 1 B- T and (1 + fi 2 ) n+1 
the equation above can be written as 



7 ,,2n+l I ,,2n+2 



XVV, + Wi_i = z,- w s 



2n+l, 



0,1, 



,2n+l, V-! 



0. 



Thus, if we require = s M to be a solution, we see that X is uniquely de- 
termined as X = —C T , where C is the companion matrix of the polynomial 
(1 + /i 2 ) n+1 . In other words, we seek a constant matrix B such that X = —C T , 
or equivalently, —M~ T = BCB^ 1 . Thus, from the basic theory of ratio- 
nal canonical forms of linear mappings, it is then clear that we can choose 
B = [v iVv • • • N 2n+1 v], where N = —M~ T and v is any cyclic vector for the 
matrix N. Especially, we can choose v = [1 • • • 0] T . The proof is complete 
when we show that the powers of N have the following block triangular form 



a—l+i—j 
i-j 







if i < j, 



(38) 
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where i, j = 1, . . . , n+1, and (N a )ij denotes the 2x2 block in the block-position 
of N a . The formula (j38|) is true for a = 0. Thus, if we assume that it is 
true for some a > 0, it follows by induction for arbitrary a: 

(^°^ = E^(^ = E(- A ) 1+ ^( :_J J h-Ar+^ 

= f a+s_1 V-A) a+i+i_: '' = ( /a+ * _ - ? V-A) a+i+i ~ : '' 

□ 

In the s-coordinates, every analytic solution of ([T]) can be expressed as finite 
sum of *-power series 

oo 
3=0 

where a M is a constant solution. 

Theorem 11. Assume that the matrix M consists of a single real Jordan block 
(| _?_?[) . corresponding to a complex conjugate pair of eigenvalues. Then, every 
analytic solution of (fJJ) is ^-analytic. 

In detail, let (/, g) be an analytic solution of (fj]) defined by f + ig = Fq + 
Fi + ■ ■ ■ + F n , where F^ is given by and the corresponding function <pk has 
power series representation 

oo 
3=0 

Then (f,g) can be expressed as a *-power series in the following way. 

n oo 

E E CmjOm,!* * = / + •" + 9^ + \ (39) 

m—0 j—Q 

where s^ is given by ( |ff7| ) and a rn ^ are constant solutions given by 

3 m 



m J m / .\ 

<W = E + f) n - J E E(- 1 ) ?C2 " m U K 

o— n i — n n V ' 



j=0 k=0 i=0 

where 



{% ■ 2k\ 

(—1) 2 I if m is even 



i 



1 .m = l/i+2fe-l> , 

|U(— 1) 2 I . I if m is odd 
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Proof. The proof is by induction over n. For n = 0, / + ig = Y]j CQj(xo + iyo) J , 



B = 



-1 




ei ( )(-A) ei 



= I. 



Thus, x = s, and since ao lAt = 6ooo = 1, the *-power series in the left hand side 
of [[59]) reduces to 

oo oo 

c ojao^ * (s M )i = Y c Oj( x o + Vo(J>)l = f + 

3=0 j=0 

Hence, the theorem holds for n = 0. If V n>tl denotes the *-power series in the 
left hand side of (3§J). Then, 



EE' 

ra=0 j=0 



(SfJl = ^C n: ja n ^ * (S^l + V n -l^(l + /U 2 ). 

3=0 



Thus, if we assume that V n -i. ^ =/ + ••• + QH> where / + ig = Fq + F\ 
• • • + F n -\, it is enough (by induction) to prove that 



^ ' c njQ<n,fj, * ( s fi)t — /« + "" + 9n^ 2n+1 , 
3=0 

where f n + \g n = F n . By letting F„ = fn + iffn denote the function F n , 
corresponding to <p(s) — s , the proof is complete if we can prove for any iV 
that 



(40) 



where a M = a n ^. We start by confirming the case when N = 1, and thereafter 
we consider the general case N > 1. 

We begin by expressing the functions fn^ and g„^ in the s-variables. Let 
z v — z o + z\fj, + ■ ■ ■ z n \i n , then 



n / . \ I 



dfi n ^ E ^ 2 J (n-IJlfl/i"-'*' 1 
f>"-<(-i)^ n _* ) . 



i\ n! 9 s 



71 / 



Thus, if we let 



D 



1 

-1 
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then, by using lcmma fTUl and the identity AD = —DA, we get in matrix notation 

t(i) " 



Jn 

o (1) 
yn 



2" 



2 n \ 



2" 
Vn 



% n 
Vn 



1=1 
2 n - l A l D 



%n—l 
-Un-l 



1=1 



•En— I 
Un-l 



— \—A n D ... - 2 n - 1 AD 2 n l] 

on L J 



f n H>("A)' 



-2 n - 1 D(-A) 2 n l]Bs 



2n+2 



- 2"B n+1J - £ 2 i_ 1 -D (- A) n ~ i+1 B itj Sj -! 

j=l \ i=l / 



2n+2 



2n+l r 

= E 

j=0 



6j 
d, 



(-A)^e lS „ 



where 



On 



»=0 



and bj = 2n\ 



n+j-1 



If we assume that n is even, since A 2 = —I, we obtain 



-M) 



2n+l 

EM) 4 



6^ 
di 



(-Ay eiSj 



r=0 
n 



+ r 



&2r 










C?2r 


— S2r+1 



r=0 



S2r 



&2rS2r 
<i2r+lS2r+l 



&2r+l 
tfcr+l 



A ei 



When n is odd, we get in a similar way 



f (i) 



r=0 



-&2r+lS2r+l 
d2 r S2r 



We now have to confirm that a u * s M = + • • • + V 2 " +1 ■ We note that it 
is enough to prove that the coefficients at the highest power of \x coincide, since 
the other coefficients are then uniquely determined up to irrelevant constant 



28 



terms. Moreover, we consider only the case when n is even, since the case with 
odd n is completely analogous. Let = X)™=o a «(l + A* 2 )*: an d express in 
the following way. 

2n+l n n 

S» = Sr ^ = X!^ 2fe + S 2k+lV)V? k = ^2( s 2k + S 2 fe+lAt)(l + - l) fe 

r=0 k=0 k=0 



3=0 k=j 

Then, we obtain 



n k 

= ^(l + ( u 2 ) fe ^a fe _ r t r 



\i=0 r=i ^ ' ) 

\r=0 j=0 ^ ' / 



It is therefore enough, for proving the case N — 1, to confirm that the following 
expression vanishes for each r = 0, 1, . . . , n. 



,-_n V / 



i=0 

2 n 



-D-)'C)"'sg ( - i ^a) ( - i ' , (T 

Thus, the case TV = 1 is proved by the following calculation. 

£<M0(3 -£<-<)(£) 

W i=0 \ * / I if r ^ k. 
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We are now ready to prove (00]) for arbitrary N > 1. We have already 
proven the case when n = 0, so it is no restriction to assume that n > 1. Let 
/i^' and g„ denote the coefficients at the lowest and highest powers of /i in 



F„ ^ by showing 



.N-l 



o-p, * (s M )* , respectively. We will prove that /„ 

that, when considered as polynomials in xq, their coefficients at the power x 2 
coincide. Thereafter we will motivate that this implies that the functions must 
be identical. 



We consider first the function F, 
follows that 



(A f ) 



Since 2^ — xq is independent of xq, it 



5^ = (z + z\[i H z n /j, n ) N = (x + (z M - x )) N 

= Xq + NxQ~ 1 (z tJi — x ) + (lower order terms in x (l.o.t.)) 
= Nx T n 



o Zp — (N — l)x v + (l.o.t. 



Thus, since (N — 1)xq is independent of /i, we get 



p(N) = [ 



n M 



— 1 

y 



d 



n—l 



(n-l)\ <V 



-z N 



i=i 



N- 



(n-l) 



]Nx Q 



/i=0 

-i 5 



n—l 



n\z N » 



+ (l.o.t.) 



/x=0 



nlx^ -1 \Nz n 



n\x 
x N ~ 



Jtf-1 i 



Nz r , 



n-l 

-E 



i=i 



Nz n , - 



NZn-l 



2 

(JV 



(a;o-iiV2/o)J+ (l.o.t.) 
l)x ] + (1-o-t.) 



l (NF^ -n\(N- l)ar ) + (l.o.t.) 



On the other hand, from lemma [TUl it follows that, so = xo + Co, where cq is 
independent of xq, and that si, . . . S2„+i are independent of xq. Thus, 



(x + (s M - x )) w = iV<~% - (JV - l)a# + (l.o.t.), 



JV-1 



N 



and therefore, using the fact that we already proved the case when N = 1, we 
obtain 

a M * (s^ = a M * (iVx^- 1 ^ — (N — l)x£ + (l.o.t.)) 



„W-1 



* a„ — (JV — 1)-T a p ) + (l.o.t.) 



< _1 (NfW -(N- l)a x 

+ ■■■ + ti 2n+1 (Ng^ -(N- l)a n x )) + (l.o.t.) 
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Hence, 



fL N} + i~9i N) = (nfV -(N- 1)(oo + ia„)*o) + (Lo.t.) 

we conclude that the coefficients of /« ' + ig^ and F,i at the power x^ 1 
coincide. 

The last part of the proof is to motivate that fn + iffn an d Fn^ must 
coincide. From the ^-multiplication theorem^ it follows that (f n , 9^) must 
be a solution of ([1]). Therefore, /i"' +i<?n = -^b + ' ' ' +F n , where the functions 
Fk have the form (fT2)l for some choice of Moreover, f n N ^ and are 

homogeneous polynomials in the variables Xq, yo, ■ ■ ■ , y n - Since Fk = Ft will 
give a polynomial solution which is homogeneous of degree r, we can therefore 
conclude that Fq — boFfi N \ F% = biF^ N \ ■ ■ •, F n = b n Fn N \ for some constants 
6o, • • • , b n . Thus, for some polynomial c which is constant in x$, we get 

/W + i~9i N) F^ = b F^ + ■■■ + bn^l + (b n 
= cx£ + Nx^- 1 (b (F^ - io) + • • • + b n -!(F^\ - x ) 

+ (b n -l){F^ -xoj) + (Lo.t.) (41) 



Finally, the polynomials Fq 1 ^ — xq,..., Fn — xq are linearly independent since, 
for each k, F^ —xq is independent of yk+i, ■ ■ ■ , y n with a non-trivial dependence 
on j/fe. Therefore, since we know that the x^ -1 -coefficient in (|4"Tj) is zero, we 
obtain that bo = ■ ■ ■ = b n -\ = and b n = 1. Hence, we conclude that f n + 
igk — F n , which completes the proof. □ 

In order to prove that any analytic solution of {1} is also ^-analytic, one 
would need to cover the last case when the matrix M consists of several real 
Jordan blocks corresponding to the same eigenvalue. This step is technically 
involved and we have omitted it here. It is however expected that the result 
holds in whole generality, as we have proved in the complex case. We conclude 
this section by stating a hypothesis for the last remaining case of the equation 

©• 

Assume that M is a real matrix with eigenvalues ±i. It is then no restriction 
to assume that M has the form M — diag(— C^ T , —C^, ■ ■ ■ , — C^ T ), where 
each block C, is the companion matrix for the polynomial (1 + /i 2 ) ,li+1 with 
n% > ri2 > ■ ■ ■ > n m . Let s~, s\, . . . , s\ n +1 , Sq, . . . , s™ nm+1 be the corresponding 
coordinates. According to theorem [4] the equation JT]) can be extended to the 
/^-dependent equation (|23|) . 

Hypothesis 1. The general analytic solution of the equation £7p, with M — 
diag(— C]" T , — C^~ T , . . . , — C ? ~ T ), is ^-analytic. In detail, every analytic solution 
{fig) of (OP can be obtained from a *-power solution = f + Vi[i + - ■ - + gfi 2n+1 
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of (ESP given by 

m 

r=l /=(ti,t a ,...,<r)6N r 

where 

(ajl := * (M ni -" 2 (4)i 2 * • ' ' * (/* np - 1_np WJr)) • (42) 

T7ie * -power should be interpreted in a similar way as the corresponding 
* -power LS6}) in the complex case. 



6 Conclusions 

The product of two holomorphic functions is again holomorphic. In terms of 
the Cauchy-Riemann equations, 



1 

-1 



this fact can be expressed through a multiplication (bilinear operation) * on the 
solution space 

(f,g)*(f,g) = (ff-g~gJ~g + gf)- 

That any holomorphic function is analytic means that any solution (/, g) of 
the Cauchy-Riemann equations can be expressed locally as a convergent power 
series of a simple (linear in the variables x, y) solution. In a neighborhood of 
the origin for example, any solution can be written as 

oo 

(f,g) = ^2( a r,b r ) * (x,y)l, where (x,y)l = (x,y) * ■ ■ ■ * (x,y) . 

The main result of this paper is that we establish similar properties for the 
more general system 

V/ = MVg, (43) 

where M is an arbitrary n x n matrix with constant entries: 

1. Any system (|43[) admits a multiplication * on the solution space, mapping 
two solutions to a new solution in a bilinear way. 



2. The analytic solutions of l]43[) are characterized by the ^-analytic functions, 
meaning that every solution of (14"3"|) can be expressed locally through *- 
power series of simple solutions. 

Systems of the form (|4"3"1) constitute only a special case of a quite large 
family of systems of PDEs which admit ^-multiplication of solutions [5]. A 
natural problem, worth studying in the future, is therefore to settle whether the 
results for the gradient equation (|43|) can be extended to more complex systems 
of PDEs. 
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